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Abstract 



The aim of this paper is to construct six - dimensional symplectic thin-lens transport 
maps for the tracking program SIXTRACK |2|, continuing an earlier report ^ by using 
another method which consistes in applying Lie series and exponentiation as described by 
W. Grobner Q and for canonical systems by A.J. Dragt [Q]. As in Ref. [0 we firstly use 
an approximate Hamiltonian obtained by a series expansion of the square root 



Furthermore, nonlinear crossing terms due to the curvature in bending magnets are ne- 
glected. An improved Hamiltonian, excluding solenoids, is introduced in Appendix A by 
using the unexpanded square root mentioned above, but neglecting again nonlinear cross- 
ing terms in bending magnets. It is shown that the thin - lens maps remain unchanged and 
that the corrections due to the new Hamiltonian are fully absorbed into the drift spaces. 
Finally a symplectic treatment of the crossing terms appearing in bending magnets is pre- 
sented in Appendix B, taking into account only the lowest order. The equations derived 
are valid for arbitrary particle velocity, i.e. below and above transition energy and shall 
be incorporated into the tracking code SIXTRACK Q. 
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1 Introduction 



Continuing an earlier report |jT| , in this paper we show how to solve the nonlinear canonical 
equations of motion in the framework of the fully six-dimensional formalism for various kinds of 
magnets (bending magnets, quadrupoles, synchrotron magnets, skew quadrupoles, sextupoles, 
octupoles, solenoids), using an approach different from Ref. namely by applying Lie series 
and exponentiation to a kick approximation |0, 3 Q 

In addition to Ref. |I| we also study the thin -lens formalism for an improved Hamiltonian 
which is exact outside solenoids and bending magnets. It is shown that the thin -lens maps 
obtained earlier remain unchanged and that the corrections due to the new Hamiltonian only 
appear in the drift spaces. 

The equations derived are valid for arbitrary particle velocity, i.e. below and above transition 
energy and shall be incorporated into the tracking code SIXTRACK 0. 

The paper is organized as follows : 

In chapter 2 the general canonical equations of motion are derived. The thin -lens method 
using Lie series and exponentiation is described in chapter 3. Using the thin-lens approximation 
the equations of motion are solved for each element in chapter 4. The improved Hamiltonian is 
introduced in Appendix A. In addition to Ref. the influence of nonlinear "crossing terms" 
resulting from the curvature in bending magnets is investigated in Appendix B and a super- 
position of a solenoid and a quadrupole in Appendix C. Finally a summary of the results is 
presented in chapter 5. 



2 The Canonical Equations of Motion 

2.1 Notation 

The formalism and notation in this paper will be identical to that used in Ref. Thus 
we will begin by simply stating the canonical equations of motion already used in this earlier 
paper and refer the reader to the latter for details. 

2.2 The Hamiltonian in Machine Coordinates 

The Hamiltonian for orbital motion in storage rings reads as [|l| : 
n{x,p^,z,p^,a,p^;s) = - [1 + f{p^)]- [1 + ■ X + ■ z] X 

( [p^ + H-zf + ip^-H-xfY^^ 



1 



[1 + f{P.)f 



[l + ■ X + K, ■ zf g ■ {z^ - x^) - N ■ xz 

+ - ■ (x^ -?>xz^) 
6 

I ^ / 4 c 2 2 I 4\ 

+— ■ [z -<ox Z + X ) 
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+ ^ ■ 7^ T ■ 7^ ■ COS 
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(2.1a) 



^There is a vast literature on solving differential equations by Lie series. A nice treatment is given by Ref. 



with 



2,2 / moc- 



Eq 



{g, N, H, Kx, Kz, X, and /i are defined in Ref. ^ ). 
Since 

\Px + H-z\ < 1 ; 

\pz — H ■ x\ <^ 1 



the square root 



■ _ [p, + H-zf + [pz-H-x] 



[1 + f{p^)r 



1/2 



in (2.1) may be expanded in a series: 



(2.1b) 



_ [px + H-z]^ + \pz-H- x] 
[1 + fiP.)? 



1/2 



[Px + H-z]'^ + [p,-H-x] 
2 ■ [1 + f{p.)? 



+ ■ 



(2.2) 



The power at which the series is truncated defines the order of the approximation to the particle 
motion. 

In the following we will use (as in Ref. [Q) the approximation: 

1 [p., + H-z]^ + [pz-H-x]^ 



n 



+ 



2 [1 + f{p.)] 

Pa-[l + K^- x + Kz- z] - f{p„ 



l[Kl + g]. x' + \[Kl - g]. z' - N . XZ + 



A 

6 



[x^ -?,xz^) + ■ iz^ -Qx^ z'' + x^) + 



24 
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13^ 2-K-h Eq 
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. 27r 
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An improved Hamiltonian is introduced in Appendix A. 

The canonical equations corresponding to the Hamiltonian (2.3) take the form: 



ds 



y 



&H 



(2.3) 



(2.4) 



with 



= {yi, 1/2, 1/3, 1/4, 1/5, ye) 

= {X, Px, Z, Pz, (T, Pa) 



and 



or, written in components : 



( S2 0\ 

0^2 

V ^2 y 



-1 

+ 1 



(2.5) 



as 



ds 



ds 



dpx 
Px + H ■ z 
[1 + f{p.)] ' 

[p^ - if ■ x] 



(2.6a) 



H-[Kl + g]-x + N-z + K,-fip,) 



[1 + f{P.)] 
_\ 

~2 



m 



dpz 
Pz — H ■ X 

[1 + fM\ ' 

+ 



H-[Kl-g]-z + N-x + Kz-f{pa) 



+A • — ^ • (2;^ — 3 2;) ; 



(2.6b) 



(2.6c) 



(2.6d) 



— a - + 
ds 



on 



dpa 

l-[l + K,-x + Kz-z]-f'ip,) 

_1 \px + H-z]^ + [pz-H-x]^ 

2' [l + fiPa)? 

l-[l + K,-x + Kz-z]-npa) 

-l-[ixr+izr]-f'iPa); 



fip^) 



(2.6e) 
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In detail, one has : 



a) 


X2 + X,V0; 














bending magnet; 


b) 


9^0; 




= K,= 


iV = 


= X = fi = H = 


-- V = 


-- : 


quadrupolc; 


c) 






= K,= 


9 = 


X = fi = H = 


V = 


: 


skew quadrupole; 


d) 


A^O; 




= K,= 


9 = 


N = 12 = H = 


-- V = 


= : 


sextupole; 


e) 






= K,= 


9 = 


N=X=H= 


-- V = 


= : 


octupole; 


f) 


H^O; 






9 = 


N = X = fx = 


v = 


: 


solenoid; 


g) 








9 = 


N = X = fx = 


H = 


: 


cavity. 



3 Description of the Thin - Lens Method 

3.1 Thin - Lens Approximation 

The equations of motion (2.6) have the general form : 



or 



with 



and 



^Vi = (l/i, 2/2, 2/3, 2/4, 1/5, l/e; s) ; (3.1a) 
(i = 1,2, 3, 4, 5, 6) 



= ^ims) (3.1b) 



= {^1, ^2, ^3, ^4, A, A) 



-[Kl{s)+g{s)]-x + N{s)-z 



2 ^ '6 



7 



-[K,{sy-g{so)]-z + N{s)-x 

+X{s) -xz-l^-iz^-Sx^z) ; (3.2d) 



1 2 2 

- ^^^^^ • + • f{p.) - ^^^l: -J/ ■ "J • /(p.) ; (3.2e) 



2 [i+/(p.r [i+/(p.)r 



Po ^0 



n ■ — ■ a + ip 

1j 



(3.2f) 



Equation (3. 1) represents a system of differential equations the solutions of which can be 
written in the form : 

y{si) = T{sf,Si)y{si) (3.3) 

by defining a transport operator T{sf,Si) connecting the final vector y{sf) at position Sf 
with the initial vector y(sj) at position Sj. 

The aim of this chapter is now to calculate the transport map T{sf, Si) (approximately) by 
using symplectic kicks. 

We achieve that in two steps : 

In a first step we decompose the r.h.s. of (3.1) into two components: 

1? = ^D + '&L (3.4) 

gathering in all terms of i) containing the external electric and magnetic fields (expressed 
by the lens functions V, g, N, H, K^, Kz, A, ). 

— * 

As a result, the component i^d in (3.4) then corresponds to the Hamiltonian 

(to be obtained from (2.3) by neglecting all external fields) leading to the (canonical) equations 
of motion for a pure drift space : 

d d ^ , . , 

-r^ ^ +^'T-D{x,Px,z,Pz,cr,Pa) 
as opx 



■'X 

Pz 



[1 + fiPa)] ' 



(3.6a) 



-^Px = —^'HD(x,Pa:,Z,Pz,a,Pa) 



— p-r = const ; (3.6b) 



8 



as opz 



■'z 
Pz 



[1 + /(p.)] ' 

j-^Pz = --Q^nD{x,Px,z,Pz,a,Pa) 
— Pz — const ; 

d d ^ , . . 

^ +-^'^D{x,Px,Z,Pz,(T,Pa) 

ds opa 



2 



as acr 

= =^ Po- = const 

(see also eqn. (2.6) ). The solutions for a drift of length I are : 



[1 + m)] 



Px Px 1 



Pi = Pz 



Pa = Pa 



The second component -di corresponds to the Hamiltonian 

'^L = 77— Ff ^ - iPx - Z - Pz - X\ + - 



2 



[1 + /(P.)] ' 2 [l + /(p.)] 



6 24 



9 



_!_ L eVjs) 
PI ' 271 ■ h ' Eo 



■ cos 



, 27r 

h ■ — ■ a + ip 

1j 



(3.8) 



containing the remaining terms in eqn. (2.3). In particular there are no or terms. 
Thus we have 



n 



and 



oy 



= 
= 

where the matrix S_ is given by eqn. (2.5). 

In the second step we replace the function -(9 in (3.4) for a thin lens of length As at position 
So by i 



with 



^D{y) + '&L{y; s) ■ As ■ 5{s - So) 
^D{y) + ^L{y; so) • As ■ 6{s - so) 



-^dy; So) 



(3.9) 



(3.10) 



and 



'HL^m So) 

H{so) 



[p,^ ■ z-p^-x] + 



H\so) 



[l + f{p.)] " ' 2 [l + /(p.)] 

[K^{so) ■ X + K,{so) ■ z] ■ f{p^) + 



2 I 2 

X + z 



\ [Kliso) + 9{so)] -x' + l [Kl{so) - g{so)] ■ z' - N{so) ■xz + 
^^''^ (x'-3xz') + ^-(z'-6x'z' + x') + 



6 

1_ _J_ eVjso) 
(3f 27T-h' Eo 



24 



cos 



■ 2tt 

h ■ — ■ a + if 
Ij 



whereby the new function ^mod in (3.9) results from the modified Hamiltonian 

Ti-mod = Hd + T-Cl ■ As ■ S{s - So) . 



(3.11) 



(3.12) 



In order to solve eqn. (3.1) using the modified function ^mod in (3.9), we then have to 
decompose the region 

As As 
■So ;^<s<so + — 



10 



of the lens into three parts : 



region I 
region II 
region III 



So 
So 



As 

^ < s < So 



€ < s < so + € ; 

As 

so + e < s < so + — 



(3.13a) 
(3.13b) 
(3.13c) 



(0 < e ^ 0) . 

As 

For region I and III we obtain a drift space of length / = — , described by the differential 

2 

equation 

^Jis) = (3.14) 
the solution of which is given by eqn. (3.7) and may be expressed by a transport operator 

The equation of motion for the central region II reads as : 

^y{s) = 4L{mso)-As-5{s-s^) (3.15) 

as 



with 



and 



(^l)^ = (^Ll, ^L2, ^L3, ^L4, ^L5, ^L6) 



d 

'&Li{y;so) = +-K—'HL{x,Px,z,p^,a,Pa) 



d 

1^L2{y;So) = --Q^T-CL{x,Px,Z,p^,a,Pa) 



, Pz Tji ^ H{sq) -X 

+X,(so) • f{pa) - [Kl{s^) + g{s^)\ ■ X + 7V(so) • z 

'^(■^o) / 2 Ksq) 3 „_^2 



(x^ - z^) - -ix^ -Zx z^) ; (3.16b) 

6 



d 

^Lz{y\so) = +^'HL{x,Px,z,Pz,a,Pa) 
dpz 

H{so) ■ X 



d - 

■dLAimSo) = -—nL{x,Px,Z,p^,a,Pa) 



(3.16c) 



11 



+K,iso) ■ f{p.) - [Kl{s,) - g{s,)] ■ z + N{s^) ■ x 

+A(so)-xz-^^-(z3-3a;'z) ; (3.16d) 
6 



d - 

^L5{y]So) = +—nL{x,p^,z,p,,a,Pa) 

OPa 



■[K^iso)-x + K,{so)-z]- f'{p^] 

1 H{sof ■ [x^ + z^] , , H{so) -Ip^-z-p.-x] 



f\p.) - tiy^}lL^^^^L^^ . f'^p^) (3.16e) 



d 

^L6(^;so) = --g^nL{x,p^,z,p^,a,p^] 



1 eV{so) . 
Po ^0 



h ■ — ■ a + if 

1j 



(3.16f) 



(see eqns. (3.10) and (3.11) ) and determines the transport map 

Tl = T(so + 0,So-0) (3.17) 

of region II. 

Finally the transport map of the whole lens takes the form : 

T(so + As/2, So - As/2) = T^(As/2) ■ ■ rB(As/2) (3.18) 

corresponding to the decomposition of the length As into three parts (see eqn. (3.13) ). 

Note that the (nonlinear) transport maps and T_i^ corresponding to (3.14) and (3.15) 
and thus also T(so + As/2, sq — As/2) in (3.18) are automatically symplectic for an arbitrary 
As due to the canonical structure of the equations of motion (see also Ref. |Q and Appendix 
AinRef. §). 

In the limit 

As — ^ 

one obtains the exact solution of the canonical equations of motion corresponding to the starting 
Hamiltonian (2.3). 

Since 2Zd is already known from eqn. (3.7) we are left with the problem of calculating the 
transport map by solving eqn. (3.15). This is done in the next section, using a Lie series 
and exponentiation. 



3.2 Integration by Lie - Series 

3.2.1 General Autonomous Case 

In the thin -lens approximation the equations to be solved are not autonomous but the 
s - dependence is trivial which reduces the calculation to an autonomous system. 



12 



An autonomous system of differential equations of the form : 



d ~ d ~ 

j-^Vi = A iVu y2, ■ ■ ■ , Vn) ; = 0; (3.19) 



= 1, 2, . . . , n) 



(no explicit s dependence) where the terms i^i{yi, 1/2, ■ ■ ■ , Vn) represent analytical functions, 
can be solved by Lie-series : 



yi{s) = e[(^-^°)^ly, (3.20a) 



with 



d ~ d ~ d 

oyi dy2 dy, 



and 



(3.20b) 



yi{so) = Vi ■ (3.20c) 



Applying this result to the canonical equations of motion : 

d 

d 

^2 = --^n{x,p,,,z,p^,a,pa) ; 
d 

^3 = +^H{x,p^,z,p,,a,p„) ; 
dp, 

d 

-^i = --Q^'H{x,p^,z,p,,a,p„) ; 
d 

^5 = +T:—H{x,p^,z,pz,(r,pa) ; 

dpa 

d 

1^6 = --^'H{x,p^,z,p,,a,p^) 



we obtain : 



l/i(s) = x{s) = e^'-'^^^x ; (3.21a) 

y2{s) = p,{s) = e(^-^°)^p, ; (3.21b) 

y,{s) = z{s) = e(^-^°)^z; (3.21c) 

y,{s) = p,{s) = e(^-^°)^p, ; (3.21d) 

y,{s) = a(s) = e(^-^o)^a; (3.21e) 

y,{s) = p^{s) = e(^-^")^p. (3.21f) 
13 



with 



and 



D 





OPx 


d 
dx 


ox 


d 

dpx 


+ 




d 

dz 


oz 


d 

dpz 


+ 


r) 


o 
oa 


' r) 
oa 


o 

OPa 


X 


= x{so) 


; Px = 


Px{so) ; 




z 


= z{so) 


; Pz = 


Pz{so) ; 




a 


= o"(so) 


; Pa = 


Pai-So) . 





(3.22) 



(3.23) 



Remarks: 



1) Using the notation of Ref. |Q, eqn. (3.20a) may also be written in the form: 

yi{s) 



if the autonomous equations of motion result from an Hamiltonian 7i. So when the approach in 
Ref. is restricted to canonical systems it is identical to the Lie Algebra method introduced 
by Dragt. 

2) Since the equations of motion (3.6) for a drift space represent an autonomous system of 
differential equations, eqns. (3.20a, b, c) can be used to determine the transport map T^, of a 
drift space. 

In this case we get by comparing (3.19) with (3.6) : 



^1 

^5 



y2 



+ ' 

; 

[1 + five)] ' 

0; 

1 - five) 
. 



^ [,.2 I ,.21 five) . 

2 ^''^''^ [1 + W 



This leads to 



Viiso + l) = e'^'vi 
14 



(3.24a) 
(3.24b) 
(3.24c) 
(3.24d) 

(3.24e) 
(3.24f) 

(3.25a) 



with 



and 



D = + + (3.25b) 



Vi = Viiso) . (3.25c) 



We then have : 



and 



Dy2 = 0; (3.26b) 

Dy^ = ; (3.26d) 

2 [1 + /(ye)] 

Di/e = (3.26f) 



D^y = for z/ > 1 



Thus: 

y{so + l) = [l+l-D]ij. (3.27) 
Putting (3.26) into (3.27), we regain eqn. (3.7). 

3) The method for calculating thin -lens transport maps described in this paper works also 
in the presence of nonsymplectic terms resulting for instance from radiation damping One 
must simply include these terms in D before expanding exp D . 

3.2.2 Calculation of the Thin -Lens Transport Map for the Central Region. 

In order to determine the transport map T^i for the central region, we investigate the 
special case : 

&,{y) = 6{s~So)-F,,{y} ; ^ = . (3.28) 



15 



Replacing the 5- function 6{s — Sq) in (3.28) by a step function of height (l/2e) and length 
(2e), we obtain in this case from eqn. (3.20) : 



with 



and 



y(s) = jexp 

for {sq — e < s < Sq + e 

y = y{so-e) 



dyi dy2 dye 



{[so + e]-[so-e])-—D 



In particular by putting s — sq + e we have : 

^(■so + e) = jexp 

= [exp[D]}y, 
which leads in the limit e — > to : 

y{so + 0) = {exp[D]}^ 



with 



y = y{so - 0) . 

Then by choosing the functions Fi{y) appearing in (3.29a) as: 

Fi{y) = i?Li(y;so) • As , 



(3.29a) 



(3.29b) 



with 'dbiiy^so) given by (3.16), one just gets the transport map corresponding to eqn. 
(3.15) in the form : 



= expp] 

as may be seen by comparing (3.28) with the r.h.s. of (3.15). 



(3.29c) 



Remark: 

The relation (3.29c) for can also be derived by solving the differential equation : 

= 4(y;so) = ^^-m , (3.30) 

which does not contain the 5 -function S{s — Sq) anymore. Writing the solution of (3.30) in 
the form 



y{s) = T{s, So) y{so) , 
16 



(3.31) 



one then obtains : 



Tl = T(so + As, So) , 



(3.32) 



as can be verified by compararing (3.29) with (3.20) and using (3.19). 

By (3.30) we see in fact that for the central region II the problem reduces to an autonomous 
one. 

For an example see Remark 2) at the end of Appendix C. 

4 Thin - Lens Approximation for Various Kinds of Mag- 
nets and for Cavities 

In this section the thin -lens transport map corresponding to the central region II (see eqn. 
(3.13) ) is calculated for cavities and for various kinds of magnets. 

4.1 Quadrupole 
4.1.1 Exponentiation 

For a quadrupole we have : 



= K, = N = X = i^ = H = V = 0. 
Then we obtain from (3.16) and (3.29b) : 

Fi{y) = 0; 

F2{y) = -fi'(so) • As • a; ; 

Fsiy) = 0; 

FM = +g{so) ■ As ■ z; 



9 ^ 



and 



FM 



0; 



Fe{y) 



. 



Thus: 



b = F2{^)- — + FS)- 



d_ 
dy4 



(4.1) 
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and 



Dy 



I \ 




\ 








Ay 



) 



with 



A 



As 



/ 

















\ 




-9 













































+9 































V 

















J 



reads as : 



M = exp 



A 



1 + A 



The transfer matrix M defined by 

y{so + o) = My(so-O) = m 



since 



and 



(4.2) 



(4.3) 



(4.4) 



DA = AD =^ D"^ = a" ^ =^ |exp z) | ^ = |exp i } ^ 



A 



for z/ > 1 . 



4.1.2 Thin - Lens Transport Map 

From (4.4) we obtain [|: 



f i 
= X 



pi = Px- 9{so) ■ As ■ x' 



^See also section A. 2. 2 in Appendix A, where a superposition of quadrupoles, skew quadrupoles, bending 
magnets, sextupoles and octupoles is investigated. 
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^ Pi + 9{so) -As-z'; 



pI = pi 



with 



y' = y(so - 0) ; 
y^ = y{so + 0) ; 



X, px, z, p^, a, Pa) . 



4.2 Skew Quadrupole 
4.2.1 Exponentiation 

For a skew quadrupole we have : 

N 



Thus we get from (3.16) and (3.29b) : 

Fi{y) = 0; 

F2{y) = +iV(5o) • As • z ; 

Fz{y) = 0; 

Fi{y) = +A^(so) • As • x; 



and 



FM 



0; 



FM 



. 



Thus: 



b = F20)- — + FS)- 



d_ 
dy4 
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and 



Dy 



( \ 

F2{S)) 








ky 



with 



As 




















\ 











































iV 





































I 

















J 



The transfer matrix reads as : 



M — exp 



A 



since 



and 



l^A 



DA = AD 



A = for i/ > 1 . 



4.2.2 Thin - Lens Transport Map 

Prom (4.8) we obtain : 



pi = pi + N{so)-As-z'; 



pf ^ pi + N{so) -As-x' ; 



pI = pi 
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4.3 Bending Magnet 
4.3.1 Exponentiation 

For a bending magnet we have : 

+ ^ 0; K,-K, = 

and 

g = Ar = A = // = i/ = y = 0. 
From (3.16) and (3.29b) we thus obtain: 



i^i( 



0; 



-[K,{so)f-As-x + K,{so)-As-f{p,) ; 



F^{y) = -[K,{so)f-/^s-z + K,{so)-^s-f{p,)- 



Thus: 



and 



with 



Fe{y) = 



D = F2(S-|^ + ^4(^)-|^ + F5(^)-|- 

5^/2 5z/4 dy5 



/ \ 



V y 



A = As- 



Dy 



Ay 



X 





















-K,-f'{p„) 




\ 

K,-f{p,) 


K,-f{p„) 


/ 
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The transfer matrix reads as : 



M — exp 



A 



1 + A 



(4.12) 



since 



and 



DA = AD 



for 1/ > 1 . 



4.3.2 Thin -Lens Transport Map 

Prom (4.12) we obtain: 

pi = p:-[K,{so)f-As-z' + K,{so)-As-f{p;); 
af ^ a'- [K, -x' + K,- z'] ■ As ■ f\pl) ; 

P/ = Pa ■ 



4.4 Sextupole 
4.4.1 Exponentiation 

For a sextupole we have : 

A 7^ 

and 

= = g = N = ^ = H = V 
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Prom (3.16) and (3.29b) we then get : 

FM = 0; 



Thus: 



and 



F2{y) ^ 

Fz{y) = 0; 



1 



--A(so)-As 



2 2 
X — Z 



FM = 0; 
Feiy) = 0. 



D = fS)-^ + fS)-^ 
dy2 dy4 



Dy 



D^y 



I \ 



FM 


V } 

=^ {exp [£>] } y 



/ \ 





V / 



= y + Dy . 



for 



4.4.2 Thin - Lens Transport Map 

From (4.15) we obtain: 



x^ = x" \ 



pi ^ pi- lx{so)-As-[{xr-{zr] 



zf ^ z' ; 



pI = pi+ A(so) -/^s-x'z' ■ 



= a' ■ 



pI = Pa ■ 
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4.5 Octupole 
4.5.1 Exponentiation 

For an octupole we have : 



II + 



and 



= = g = N = \ 
Then we obtain from (3.16) and (3.29b) : 

FM = 0; 



H = V = 



--/x(so) • As 



0; 



--/x(so) • As 



— 3 xz^ 



0; 



Thus: 



and 



D = F2(|)- — + F4(^)- — 
dy2 dy4 



Dy 



( \ 

F2{y) 







V ) 



( \ 






V j 



for V > \ 



(4.16) 



(4.17) 



{exp [-D] } y = y + Dy . 



(4.18) 
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4.5.2 Thin -Lens Transport Map 

Prom (4.18) we obtain: 

= x" ] 



pi 



pi - -i^{s^)-^s-[{x^f-2>{x') {z^)\ 



pI = - g/^(so) • As 



'Jf - 3 ix'? f/l 



pI = pI 



4.6 Synchrotron - Magnet 
4.6.1 Exponentiation 

For a synchrotron magnet we have : 

^ ^ 0; Kl^Kl ^ with K,-K, = 

and 

AT = A = // = ii" = 1/ = 0. 
We thus obtain from (3.16) and (3.29b) : 

FM = 0; 

F2{y) = -G,{so)-^s-x + K,{so)-^s-f{p,) 
FM = 0; 

F,{y) = [K,{so)-x + K,{so)-z]-As-f{p,); 
FM = 0; 

(Gi = Kl + g- G2 = Kl-g) 
25 



Thus: 



and 



with 



dy2 dyi dy^ 



( \ 


FS) 

fS) 

\ ) 



Dy = 



= Ay 



As 



V 

The transfer matrix reads as : 























\ 


-Gi 
















































— G2 










• fiPa) 









■ f'iPa) 



































/ 



M — exp 



A 



since 



and 



= 1 + A 



DA = AD 



A = for 1/ > 1 . 



4.6.2 Thin - Lens Transport Map 

From (4.23) we obtain : 

pi = pi-Gi{so)-As-x' + K,{so)-As-f{pi,); 

pi = pl-G2{so)-As-z' + K,{so)-As-f{pl)- 
af = a'-[K,{so-x + K,{so)-z]-As-f\pl)- 

pi = Pa ■ 
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4.7 Solenoid 
4.7.1 Exponentiation 

For a solenoid we have : 



H ^ 



and 



= K, ^g^N = X = ii = V = 0. 
Using then (3.16) and (3.29b) we obtain: 



Fi{y) 



F,{ 
Fei 



, H(so) ■ As 
+ Tr—r, — ^ • z ; 



[1 + fiP.)] 
Hjso) ■ As 

Hjso) ■ As 

Hjso) ■ As 
■[1 + /(P.)] 
H{so) ■ As fip^ 



■ \Pz - H{so) ■ x] ; 



X 



■ \Px + H{so) ■ z] ; 



[l + /(p,)] [l + /(p.)] 12 



- H{so) -[x +z] + \p^- z-p^ 



. 



Thus: 



D = fS)-^ + fS)~ + fS)~ + fS)-^ + fS) 
dyi dy2 dy^ dy^ 



and 



D 



f X \ 

Px 
Z 

V Pz J 



( F,{^) \ 
F2i3) 

fS) 
\ fS) J 



A, 



( X \ 

Px 

Z 

\ P J 



Da 



Dpa = =^ |exp D ] Pa = Pa 
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with 



As- 



[1 + fiPa)] 



/ +H \ 

-H"^ +H 

-H 

\ -H -H^ / 



We decompose the matrix Aq into the components Aqi A02 '■ 



with 



As 



[1 + m)] 



and 



A02 = As- 



H 



[1 + fiPa)] 



/ 








\ 




-1 



















V 





-1 


0/ 








+1 














+1 


-1 














-1 









The transfer matrix for 



yo 



f X \ 

Px 

Z 

\ P J 



reads as : 



Mo = exp A 



smce 



DAo = AqD =^ D'Vo = Alo ^ {exp[i)]}^o = {exp[io]} 



Using the relations : 

Aoi ■ A02 — A02 ■ Aoi 



exp 



exp 



A 



01 



- exp 



A 



02 



and 



for i/ > 1 



exp 



A 



01 



1 + A 



01 
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as well as 



we get : 



/ 








1 


\ 













1 




-1 











V 





-1 





0^ 


/ 








1 


\ 













1 




-1 











V 





-1 





0^ 



2n 



- (-ir-1; 



2n+l 



exp 



A 



£102 











1 


\ 













1 




-1 











v 





-1 





0^ 



with 



Therefore : 



y ^ 

nito (2n)! 



(-1)" ■ (A0) 



2n 



+ E 



% (2n + l)! 



-ir-(A0) 



2n+l 



/ 








+1 





\ 













+1 






-1 













V 





-1 








/ 



= l-cos(Ae) + 

/ cos(Ae) 



/ +1 \ 

+1 

-10 

V -1 / 



sin(Ae) 



V 



+sin(Ae) \ 

cos(Ae) +sin(Ae) 

-sm(Ae) cos(Ae) 

-sin(Ae) cos(Ae) / 



AG 



//•As 

[1 + f{Pa)\ 



01 



/ cos(Ae) +sin(Ae) \ 

cos(A0) +sin(Ae) 

-sin(Ae) cos(Ae) 

V -sin(Ae) cos(Ae) / 



For the variable a we get from (4.26b) : 
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r-, , ..^ ^1 {ys ■ - [H ■ yi - Va] ■ - yi ■ - [H ■ y:i + y2\ ■ — 



f (-//) • f{p,) ■ As 

I [1 + m)f 

H-As {-H).f'ip,).As 



1 



H -0)1+ ys) + (^2 -ys-m- yi) 



[1 + /(p.)] [1 + f{Pa)f 

X {^3 • [-ff • yi - iji] -[H -yi- y^] -ys-yi-iH ■y3 + h] + [H -113 + y2] ■ yi} 



= 0; 



(4.34) 



Remarks: 



1) We have used a separation of the matrix into the components Aq^ and 74o2 with 



and 



^01 ■ =^io2 — Ao2 ' Aoi 



Af)i nilpotent . 



Since in addition Aq2 is diagonaUzable, eqn. (4.28) represents an additive Jordan decomposi- 
tion. 

2) The decomposition (4.28) factors Mq into a rotation map and a focussing map. 



4.7.2 Thin -Lens Transport Map 

Equations (4.31), (4.34) and (4.26c) finally lead to: 

yj = Moyo'; 



(4.35a) 



a' — a 



f 



Pi 



with 



• Ae 



[1 + !{Pl 



+ 



Ae 



H{so) ■ As 
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(4.35b) 
(4.35c) 

(4.35d) 



and Mq given by (4.33), where the matrix Aqi appearing in (4.33) takes the form: 

ioi = AQ-H{so) 



/ \ 

-10 



V -1 y 



(see eqn. (4.29a) ). 

In Appendix B the superposition of a solenoid with a quadrupole is investigated. 



4.8 Cavity 

4.8.1 Exponentiation 

For a cavity we have : 



V ^ 



and 



^ ^ g ^ N ^ X 
Then we obtain from (3.16) and (3.29b) : 
Fi{y) = 0; 



^ l_i ^ H ^ 



Thus: 



and 



Fei 



0; 
0; 
0; 

1 eV(so) 



sm 



D ^ F. 



n ■ — ■ a + (f 



d 



As 



dye 



Dy 









V Fe{^) J 



D^y 



( \ 





V / 



for V > \ 
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{exp[b]}jj = ij + bj). (4.38) 



4.8.2 Thin -Lens Transport Map 



From (4.37) and (4.38) we obtain : 



f i 
X-' = X 



pi = P1 



zf = z' 



pi = pI ; 



f 1 eVjso) . 

Pi = P„ + -^-—^-sm 



\ 27r . 

h - — ■ a' + 1^ 

1j 



As 



5 Summary 

As a continuation of Ref. |[|, we have shown how to solve the nonlinear canonical equations 
of motion in the framework of the fully six-dimensional formalism using Lie series and exponen- 
tiation. Various kinds of magnets (bending magnets, quadrupoles, synchrotron magnets, skew 
quadrupoles, sextupoles, octupoles, solenoids) and cavities are treated, taking into account the 
energy dependence of the focusing strength. 

In Appendix A we have introduced an improved Hamiltonian (Appendix A), which is exact 
outside bending magnets and solenoids. 

Since the equations of motion are canonical, the transport maps obtained are automatically 
symplectic. 

The equations derived are valid for arbitrary particle velocity, i.e. below and above transition 
energy, and have been incorporated into the computer program SIXTRACK. 



Acknowledgments 

We wish to thank D. P. Barber for very stimulating and interesting discussions and for 
carefully reading the manuscript. 



32 



Appendix A: Thin - Lens Approximation with an Improved 
Hamiltonian 



In chapter 4 we have used an approximate Hamiltonian which is obtained by a series 
expansion of the square root 



The aim of this Appendix is to repeat the calculations of chapter 4 in the absence of solenoids 
with an improved Hamiltonian by using the uncxpanded square root. The new Hamiltonian 
is again decomposed into a lens part Ti^ and a drift component Tio- Then we apply the 
thin -lens approximation described in chapter 3 by defining the modified Hamiltonian 



(see eqn. (3.12) ) and solve the equations of motion resulting from the improved Hamiltonian 
for the regions I, H, and HI of eqn. (3.13). It is shown that the thin- lens maps for the central 
region II obtained earlier remain valid and that corrections induced by the new Hamiltonian 
only appear in the solutions of the drift spaces (regions I and HI). 

The improved Hamiltonian introduced in this Appendix is exact outside the bending mag- 
nets. Inside a bending magnet we neglect nonlinear crossing terms resulting from the curvature. 
These terms are investigated in Appendix B. 

A.l The Improved Hamiltonian 

In the absence of solenoids the Hamiltonian for orbital motion in storage rings takes the 
form : 




up to first order in terms of the quantity 



Hmod = Hd+Hl- ^S- 5{s - So) 




+--[l + K^-x + K,-zf---g-{z^ 
+^ ■ {x^ -Sxz'^) 



I ^ / 4 c 2 2 I 4\ 
H 7 ■ (Z —OX Z + X ) 



x^) — N ■ xz 




1 L eV(s) r, 27r 




(A.l) 



(see eqn. (2.1) ). 
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If we neglect nonlinear crossing terms containing the factors Kr^ ■ x and K^- zwe may write : 



2 , 2 ^ 1/2 



[1 + /(p.) 



2 2 ~) 1/2 



Then we obtain : 



[1 + /(P.)]^ 



2 , 2 ^1/2 



^ [1 + f{Pa)? . 

[K,-x + K,-z\-f{p,) 



+ \[Kl + g]-x' + \[Kl-g].z'-N.xz 



A 
6 



+ ^ • (x^ - 3 xz"") + ^.(z^-6x''z^ + x^) 
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^ L eV{s) 



n ■ — ■ a + (f 

Ij 



(constant terms which have no influence on the motion have been dropped). 
In particular we get for the Hamiltonian of a drift space : 



'HD{x,Px,z,p^,a,p^;s) = Pa - [1 + /(Pa)] ■ i 1 



2 , 2 ^1/2 
Px+Pz \ 



[1 + f{P.]' 



and the lens component of Ti : 



(A.2) 



(A.3) 



D 



then reads : 



nL{x,Px,z,p^,a,p^;s) = -[Kx- x + K^- z]- f{pa) 



+ l[Kl + g]-x' + \[K'^-g]-z'-N.xz 



X 



3 o^„2^ , V_6^2^2^^4) 



+ -.ix^-3xzn + ^^ 



+ 



1 L eV{s) 



f3^ 271 -h Eo 



cos 



\ 27r 

n ■ — ■ a + (p 

Ij 



(A.4) 
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A. 2 Equations of Motion 
A.2.1 Drift Space 

The Hamiltonian (A. 3) for a drift space leads to the canonical equations of motion: 

_d _ dHp 
ds dpx 



2 , 2 ^ -1/2 
Pi+Pz \ Pec 



d_ _ dUp 
ds dx 



— Px — const ; (A. 5b) 



ds dpz 



2 , 2 ■\ -1/2 
Pi+Pz 1 Pz 



Pz = const ; (A.5d) 



ds 



a 



dHp 
dpa 



2 , 2 ^1/2 



[1 + fiP^W 



ri ^ ffr. ^^ 1 /i Px. + p^ r'^' npl + Pl] ^ 



9 I 9 
P7r + P. 

+ 



1/2 



pI+pI 



'1/2 



[Px + pI] 

[1 + /(p.)]^ 



/'(p.) 
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9,9 ^ -1/2 



(A.5e) 



Pa 



dHp 
' da 



Pa = const 



(A.5f) 



The solution of eqn. (A. 5) for a drift space of length/ reads as : 



= x' + <1 



Px Px ' 



= Z' + {1 



f 

Pz 



a' 



f 

Pi 



= Pz 



pi 



[1 + /(p^)]^ 



[1 + /(p^)]^ 



[i+/(p^)] 



-1/2 



[1 + /(P^) 



1 - /'(p:) 



• Mf^{plf\ \ 

. [i+/(p^F J 



-1/2- 



= p; 



(A.6a) 
(A.6b) 

(A.6c) 
(A.6d) 

(A.6e) 
(A.6f) 



Remarks: 1) Using (A. 5a) and (A. 5c), eqn. (A.5e) may also be written in the form: 

d 



ds 



<y = l-f'{Pa)-Jl + ix'y + izr 



This result can be obtained directly from the defining equation for a : 

s -vo- t{s) ; 



with 



and leads to 



a 



d ^ dt 

ds ° ds 



dt = - • Vds^ + dx^ + dz^ 

V 



(A.7f) 



(A.8f) 
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which agrees with eqn. (A. 6), since 



(see Appendix C in Ref. ). Since 



f {Pa) = — 
V 



1 



one may write : 

This approximation was used in Ref. 

2) As in Ref. [|I|] one obtains for a drift space : 

x'{s) = const. =^ x{s) = x{so) + x'(so) • (s — sq) ; 
z'{s) = const. =^ z{s) = z{so) + z'{so) ■ (s — sq) , 

(see eqns.(A.5a, b) and (A.5c,d) ) but the connection betwenn y' and Py {y = x, z) is modified 
(see eqns.(A.5a, c) and (3.7a, c) ). 

A. 2. 2 The Central Part 

The equation of motion for the central part (region II in eqn. (3.13) ) due to the Hamiltonian 
(A. 2) reads as : 

^y{s) = F{y)-6{s-So) (A.9a) 
as 



with 



and with ^L{y', s) given by: 



F{y) = i3L{y;so)-As (A.9b) 



d 

^Li{y;s) = +—nL{x,P:„z,p^,a,Pa;s) 

OPx 



= ; (A. 10a) 

d 

^L2{y;s) = -—nL{x,p^,z,p^,a,p„]s) 

= +KM ■ f{P^) - [Kl{s) + g{s)] ■ X + N{s) ■ z 
2 6 



d 

^L3{y]s) = +—nL{x,p^,z,p^,a,p^;s) 

OPz 

37 



= 0; 



(A.lOc) 



d_ 

dz 



= +K,is) ■ fip^) - [Klis) - gis)] ■ z + N{s) ■ x 



+A(s) • xz 



3x^z) 



(A.lOd) 



d 

i^L5{y;s) = +^—'HL{x,Px,z,p^,a,Pa;s) 

OPa 

= -[K,{s)-x + K,{s)-z]-f{p,) ; 



(A.lOe) 



T<l{x,Px, z,p^, a,p^; s) 



1 eV{s) 



sm 



Eo 

resulting from the Hamiltonian Ti^ in eqn. (A. 4). 



n ■ — ■ a + (p 

1j 



(A.lOf) 



Since the relations (A.lOa-f) coincide with eqns. (3.16a-f) in the absence of solenoids 
[H = 0), the thin- lens transport maps calculated in sections 4.1-4.7 remain valid also for the 
Hamiltonian (A. 2). Thus the corrections resulting from the new Hamiltonian are fully absorbed 
in the solutions for the drift space as can be seen by comparing (A. 6) with (3.7). 

As an example we consider the superposition of quadrupoles, skew quadrupoles, bending 
magnets, sextupoles and octupoles and obtain from (A. 10) : 







Kxiso) ■ fiPa) - Gi{So) ■ X + iV(so) • z 



■^■{x'-3xz')]-As; 



0; 



G2{sq) ■ z + iV(so) • X + A(so) ■ xz 



6 



z-" -3x^z) 



As; 



[Kx{so) ■ X + K,{so) ■ z] ■ f{pa) ■ As ; 



0; 
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Thus: 



O = + (A.11) 



We then have : 



L>yi = 0; (A.12a) 

Dh = [K,{so)-f{Pa)-G,{s,)-x + N{so)-z-^-{x^-z^) 

{x^ -?,xz^)\ ■ As ■ (A. 12b) 



6 

L>y3 = 0; (A.12c) 



Dy4 ^ {Kz{so) • fipa) - G2{so) ■ z + N{so) ■ x + A(so) • xz 

{z^ -2,x^ z)\ ■ /^s ] (A.12d) 



and 



Thus: 



6 

by^ = [K,{so) ■ X + K,{so) ■ z] ■ f\p,) ■ As ; (A.12e) 
^^6 = (A.12f) 

b^y = for i/ > 1 . 

^(^0 + 0) = ^ + b^ (A.13) 



with 



— * 

y = y{so) 



and by given by (A. 12). 



Equation (A.13) contains as special cases the transport maps of simple quadrupoles, skew 
quadrupolcs, bending magnets, sextupoles and octupoles which are identical with those already 
derived in section 4. 



Remarks: 



1) As in chapter 4 the transport maps 

y [so-^ As^ — > y (so + ^ As^ 
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described by a composition of (A. 6) and (A. 13) (combining the regions I, II and III in eqn. 
(3.13) ) are symplectic for an arbitrary As due to the canonical structure of the equations of 
motion. Furthermore one obtains the exact solution corresponding to the Hamiltonian (A. 2) 
for As — > . 

2) The Hamiltonian (A. 2) is exact for a straight section with 

i.e. outside the bending magnets. For a bending magnet with 

+ + ; K^-K, = ^ 

the exact Hamiltonian reads as : 

Hftend = p.-[l + /(p.)]-[l + i^.-X + i^,-^]- 1 ^ 



[i + Zb.)]^ 



1 r .9 1 



nn-[l + fiPa)] • [X. • X + X, • ^] • |l - + /(p V 1 
-f(p,) ■ [if. • X + • 2] + i if,^. + i ■ 2^ (A.14) 



with 



{2 I 2 1 1/2 

The drift part Hd of the bending Hamiltonian coincides with eqn. (A. 3). Thus for the 
regions I) and HI) we have again the solution (A. 6). 

For the central part described by the equations : 

d - 

^Ll{y;So) = +-K—'^bend{x,Pa:,Z,Pz,a,p„) ] (A. 16a) 

d ^ 

^L2{y;so) = --Q^Hbend{x,Pa^,z,p^,a,p^) ] (A.16b) 
d 

^L3{y;So) = +^'Hbend{x,Px,Z,Pz,a,Pa) ] (A.16c) 
OPz 
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d 

^LA{y]So) = --^'Hbend{x,px,z,pz,a,pa) ] (A.16d) 
d 

^L5{y;so) = +-—nbend{x,p^,z,p,,a,p^); (A.16e) 
d 

^L6{y;so) = --^'Hbend{x,px,z,p^,a,pa) (A.16f) 



(see eqns. (3.10) and (3.11) ) we obtain the Hamiltonian : 

'Hbend = -[1 + f{Pa)] ■ [Kcc{So) ■ X + K^{so) ■ Z] 



2 , 2 ^1/2 
Px+Pz 1 



[1 + fip^)r 



-fiPa) ■ [K^so) ■ X + K,(so) ■z] + ]^ [K,{s,)f ■ + i [K.^{s,)f ■ . (A.17) 

Using (A. 9b) and (3.29), we can calculate y{so + 0) up to an arbitrary order by a series 
expansion of exp [D] . This shall be done in Appendix B where it is shown how to express the 
transport map by elementary functions (without disturbing the symplecticity), restricting the 
crossing terms to the lowest order and using a special decomposition of Idhend into two parts. 



Appendix B: Bending Magnet with Nonlinear Crossing 
Terms Resulting from the Curvature 

B.l The Hamiltonian 

By expanding the Hamiltonian (A. 14) and keeping the nonlinear crossing terms of lowest 
order resulting from the curvature , the whole Hamiltonian of a horizontal bending magnet 
{K^ 7^ 0; Kz = 0) can be written in the form : 

'Hbend = 'Hd+'Hl (B-1) 

with a lens component 

= -K^{s) ■ X ■ f{p^) + \ [K^{s)f -x' + Y ^j'l^^^ ■x-[pI + P% (B.2) 

and a drift component I-Ld given by (3.5) or more precisely by (A. 3). 

In order to obtain a thin -lens approximation of a horizontal bending magnet due to the 
Hamiltonian (B.l) which can be expressed by elementary functions, we decompose the lens 
component Tii, of eqn. (B.2) into two parts: 



with 



Ul = nf + nf (B.3) 



= -K.{s) ■ X ■ f{p^) + \ [K^{s)f ■^' + \- ^^ ^j'l^^^ -x-pI; (B.4a) 

(2) 1 K,{s) 

2 [1 + f{p^)] • ^""-^"^ 



■^A vertical bending magnet can be dealt with analogously. 
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For each part we then separately construct a thin lens (placing one lens behind the other). 
This is achieved replacing TCbend in (B-l) by 

Hbend = nD + S{s-[so-0])-As-ni^ + 5{s-[so + 0])-As-n^^^ (B.5) 

with 

= -K.(.o) ■ X ■ f{p^) + \ [K^{so)f -x' + l- • ^ • (B.6a) 

- 2 [1 + /(p.)] ^ ■ ^^-^^^ 



B.2 Thin - Lens Approximation 
B.2.1 The Component f^l^ 

Prom eqns. (3.29b) and (B.6a) we have: 



d 

OPx 



= 0; 



(B.7a) 



F2{y) = 



9 



-As • ^ T^S^^^ (x, p^, z, p„ a, p^) 

^s-\^-[K,{so)f-x + K,{so)-f{pa) 
d 

+As • — n^l^ {x, p^, z, p^, a, p^) 

[l + fiPa)]' ''''''' 

-As ■ -^H^L\x,Px,z,Pz,a,Pa) 



1 K,{so) 



2 [1 + /(P.)] 



= 0; 



(B.7b) 



(B.7c) 



(B.7d) 



+As • -^Hl^ {x, Px, z, p^, a, p^) 

OPa 



2 [l + /(p.)] [l + /(p.)] 



-As ■ -^n'i\x,px,z,p^,a,pa) 



= . 



; (B.7e) 



(B.7f) 
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Thus: 



and 



We then get : 
a) For X : 



d d d 

D = F2{^)--r + Fs{^)- — + FS)- — 
dy2 dy3 dy^ 



/ \ 



fS) 
V ) 



Dy 



Dx ^ 



Thus: 



D^x = for i/ > 



{exp[D]} 



OO -j^ 

X = ^2 ~\ ' X — X 



n = 



b) For 



F>Px = 



F2{y) 

As ■ i-[K,{soW ■ X + K,{so) ■ f{p.) - 



1 K^jsp) .2 
2'[1 + /(P.)]'^^ 



As ■ -[K,{so)Y ■ m + K,{so) ■ fim) - o 



1 K,{so) 



2 [l + ZM 



D'p, = As-D -[i^,(so)]'-yi + i^.(so)-/(y6)-7^ 



1 KM 



2 [l + /(y6)] 



•2/4 



0; 



=^ D'^p^ = for i/ > 1 . 

Thus: 

{expp]} 



y2 + As-{ -[Kx{sa)Y ■ yi + X^(so) • /(ye) - 



1 KM 



Px + As-l -[KM? • ^ + KM ■ fipa) 



2 [i + f{m)] 
1 KM 

2'[1 + /(P.)] 



■2/4 



■P2 
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c) For z : 



D^z = for 1/ > 1 . 



Thus: 



{expp]} 



d) For pz 



Dpz = 0; 



Thus: 



D'^pz = for 1/ > 



{exp [D]} 



= Pz 



e) For cr 



Da = 



= As • <^ -K^{so) ■ X ■ f'{p„) - 



1 K^{so) f'iP.) 



^1 ^4 



2 [i + ZM [i + ZM 



= 0; 



=^ D'^a = for z/ > 1 
Thus: 

{exp [b]} 



a = a + D a 
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1 




five) 


~ 2 


[1 + fim)] 


[1 + fim)] 


1 




f'im) 


~ 2 


[1 + fm 




1 







yivl 



■yl 



f) For : 



Dpa = 0; 
b^'p^ = for 1/ > 



Thus: 



(B.14) 



{exp[D]}p, = . (B.15) 



Equations (B.10-15) finally lead to: 

= x' ; (B.16a) 

pi = + As • l^-[K^{so)r ■ + K^so) ■ fipl) - \ ■ • ^Vlf } ; (B.16b) 

= + A.. (B.16C) 

pI = (B.16d) 

- .<.A.....{-^.W.m,-l.^.^.„,^}; (B.iae) 

P/ = Vl- (B.16f) 



B.2.2 The Component U^^ 

Prom eqns. (3.29b) and (B.6b) we have: 

Fi{y) = +As- —n'i\x^Px,z,Pz,(T,Pa) 

OPx 
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K^jsp) ■ As 

A ■ [yi y2] ; 



_1 K^{so) ■ As 2 

2 [l + fiPa)] 

+As ■ — 'H^L\x,p^,z,Pz,a,Pa) 

OPz 



= 0; 



•9 ^7(2), 



-As • — Hl {x.po,, z,pz, a, Pa) 



0; 



+As ■ — iif^ (x, p^, 2;, p;,, (7, p^) 
1 X.(so) • As f{pa) 2 

-1a 



2 [1 + /(p.) 



■ bi 2/2] 



9 



-As ■ ^ h'l^ {x, p^, z, p„ a, Pa) 







whereby we have written for abbreviation : 



A - 



i^.(so)-As 
+ 



Thus: 



and 



oy2 dy5 



dyi dy2 2 [l + f{pa)] 



F2(y) 





fS) 
V J 



yiy2 



Dy = 
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We then obtain : 



a) For X : 



bx = = A■[y^y,] ; 



b'^x 



= bF,{y) 



D^x 



1 -^2 \ 



1 
2 



d 



d 



1 



= • ^ {[yi y2] ^2 - ^ ^2 • 2 [yi u] 

= 0; 



Thus: 



i:)'' X = for i/ > 2 . 



{exp[L>]} 



n = 

1 yl^ 

= yi + >1 • [^1 ^2] + - • — yi vl 



(B.20) 



b) For : 

A 

Dpx = F2{y) = ---yl; 
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Thus: 



D'^p, = DF2{y) 



d 



2) v^-- 



2" 

'2, 



d_ 
dm 



d 



-2) -'-Y^-wJ'"' 



A 



{exp[L>]}p^ 



00 -| 



n = 



n = 



y2 



1 + 4^2 



(B.21) 



c) For z : 



Dz ^ ; 



Thus: 



D'^z = for i/ > 



|exp [-D] I ;S = z . 



(B.22) 



d) For pz : 



Dp, = 0; 
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D^'p, = for 1/ > . 



Thus: 



{exp[L>]}p^ = Pz ■ 



(B.23) 



e) For cr: 



Da 



[yivi] ; 



1 



4 + 



2 [1 + /(p.)] 

2 ■ [1 + /(p.)] 
= 0; 



[yi ^2] • j/2 - ^ ^2 • 2 [j/i ^2] 



i/i ^2 



D'' 0- = for 1/ > 1 . 



Thus: 



{exp[L>]}a 



C» 1 



n = 



n! 



= a + Da 

1 . /'(Pa) r. .21 



XPx 



(B.24) 



f) For : 



Dp, = 



D^'p^ = for i/ > . 



Thus: 



{exp[D]} 



Pa = Pa 



(B.25) 
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Prom eqns. (B.20-25) we finally get : 



pi 



Pa 



with A given by (B.17). 



x'-\l + A-pl + -A'-{p^ 



i\2 



x'.\l + -A-p 



Pi 



Pi = p. 



Z 1 



(j'--A 



Pa 



f'iPa) 



(B.26a) 

(B.26b) 
(B.26c) 
(B.26d) 

(B.26e) 
(B.26f) 



Remarks: 



1) The thin -lens transport map resulting from eqn. (B.5) reads as: 



Tr 



exp 



exp 



^(2) 



whereby D^^^ corresponds to the component Hi^^ and D^'^^ to Hl^. One could also use the 
more symmetric composition : 



Tl = exp 



2 



• exp 



£)(2) 



• exp 



2 



2) As mentioned in the remark at the end of section 3.2.2, one can calcTilate the transport 
map alternatively by solving the differential equation (3.30). To illustrate this method we 

^ (2,) 

use the lens corresponding to TC^ . In this case the differential equations (3.30) take the form: 



X 



Px 



z' ^ 0; 



A 






■ Px ) 


1 


A . 


~2 




0; 





2 . 

x 1 



(B.27a) 

(B.27b) 
(B.27c) 



= 0; 
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(B.27d) 



(B.27e) 



Pa' ^ 0, 



(B.27f) 



which are solved by : 
x{s) - 

Px{s 



z{s 



Pz{s 



a{s 



Pais 



x{so) ■ 



1 A "'^ 

Px{so) 



1 A 



= Pz{so) ; 



/ N 1 ^ f ba(go)] ^ N 2/ N / N 



2 As l + /baM 



Pa{so) 



Using then the relations : 

y' = y{so) ; 

y ^ = y{so + As) , 

(see (3.31) and (3.32) ) one indeed regains eqns. (B.26a-f). 



(B.28) 
(B.28a) 

(B.28b) 
(B.28c) 

(B.28d) 
(B.28e) 



Appendix C: Superposition of Solenoids and Quadrupoles 

C.l Exponentiation 

For a superposition of a solenoid and a quadrupole we have : 

H ^ 0; g ^ 



and 



Prom (3.16) and (3.29b) we then obtain: 

, H{so) ■ As 

^^^^ ^ ^¥vjm-'' 

FM = + ^}'']]^' -\p.-H{so)-x\-g-x.As- 



(C.l) 
(C.la) 
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Hjso) ■ As 
Hiso) ■ As 



X 



H{so) ■ As fip,) 



\Px + H{sq) ■ z\+g- z- As ] 



- H{so) -[x + z ] + \p^ - z - p^ 



. 



Thus: 



and 



t, = f,(J).|_ + f,(s).|_ + F3te).|- + F.(S).|- + F.te) 



D 



f X \ 

Px 
\Pz J 



( Hy) \ 

F20) 
V F^) J 



A, 



( X \ 

Px 

V P ) 



with 



and 



Da = F5( 
Dpa = 



{exp [i)] } = 



= As 





^ 





+H 


1 


-[H' + g] 










-H 










\ 


-H 









5 



[1 + /(p.)] 

We decompose the matrix Aq into the components A^i and A 02 ■ 



Aq — Aqi + y4Q2 



with 



Aoi = As- 



[1 + fiPa) 



/ \ 

-10 



V -1 y 
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and 



L02 



As- 



H 



[1 + /(p.)] 



/ +1 \ 

-{g/H) +1 

-10 

V -1 +{g/H) Oy 



(C.7b) 



Since 

the transfer matrix for 



D^Vo ^ A^Vo =^ {exp[L>]}yo = {exp[io]}yo 



/ ^ \ 

z 

\ P J 



(C.8) 



reads as : 



Mo = exp 



4o 



(C.9) 



Using the equations 



we have : 



Aoi ■ Ao2 — Ao2 ■ Aoi 



Also since 



we get : 



Because 



exp 



A 



exp 



Aoi 



■ exp 



Ao2 



Aoi 



for u > 1 



exp 



Aoi 



1 + Aoi 



( 








1 


\ 




-{g/H) 








1 




-1 











\ 





-1 


Hg/H) 


I 


( 








1 


°\ 




-{g/H) 








1 




-1 











\ 





-1 


Hg/H) 


0/ 



2n 



(-ir-1; 



2n+l 








1 




-(^/^) 








1 


-1 











V 


-1 


Hg/H) 


/ 
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we obtain : 

exp Aq2 = 



(-!)"• (Ae)^"-i 



+ E 



% (2n+l)! 



1 • cos(Ae) + 



-1)" • (Ae) 



2n+l 



/ 





+1 


\ 


-i9/H) 








+1 


-1 











V 


-1 


+{9/H) 


/ 









+1 





\ 


(g/H) 








+1 


-1 
















-1 


+ i9/H) 





/ 



sin(Ae) 



/ cos(Ae) 

-(^/i/) ■sin(A0) cos(Ae) 
-sin(Ae) 

-sin(Ae) +{g/H) -smiAQ) cos(Ae) / 



+ sin(Ae) \ 

+sin(Ae) 
cos(Ae) 



with 



AO = 



Hjso) ■ As 
+ 



Then from (C.9) : 



Mo = 



X 



l + Aol 



( cos(Ae) 

-{g/H) ■ sin(Ae) cos(Ae) 

- sin( AG) 

V -sin(A0) 



For the variable a we obtain from (C.3b) : 



+ sin(Ae) 


cos(Ae) 



\ 
+ sin(Ae) 





. (C.IO) 



-{g/H) ■ sin(Ae) cos(Ae) / 



[fS) ■ 4 + fS) ■ 4 + fS) ■ 4 + fS) ■ 4) FM 

-ig/H).y,.^^+ig/H).y,.n 



-H) ■ f{p„) ■ As 
[1 + f{P.)f 



l^H ■{yl + yl) + {y2 ■ m - ■ Vi) 
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H ■ As i-H) ■ f'ip^) ■ As 



X {m -[H -yi- Va] -[H -yi- y^] ■ ys - yi ■ [H ■ y3 + ^2] + [H ■ys + ya] • yi 

-id/H) ■ [yi -m + m-yi]} 



-2 



H ■ As i-H) ■ fip,) ■ As 



■ {g/H) -yi-yz 



2 ■ [AG]' 



[1 + m.)] 

Furthermore, using the relations : 

D^^[yim] = 
D'^'^'iyiys] - 



{g/H) -yi-m 



-1)" • [2 • Ae]2"+^ • 1 [yl - yl] ; 



00 1 

{exp[D]}[y,y,] = ^ ■ D'' [m 



n = 

00 



1 



00 -j^ 



E 



E 



(2^)! 



(-ir-[2.Aep".[yiy3] 



+ E 



^^^(2;^-(-ir-[2-Ae]--.i[^i-^?] 



and 



with 



[yi ya] • cos[2 AO] + - [y| - yl] ■ sin[2 AO] 



DF,{^) = «:-£''[yiy3] 

1 /'(p.) g 



K 



2[1 + /(P.)]' ^ 

1 f'iPa) g_ 

2[l + /(p,)] if' 



we get : 



00 -j^ 

{exp[D]}a = E --Da 



n = 



n = 2 
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OO -1 



n = 2 



OO 1 



n = 2 



OO 



n = 2 ^! 



OO 2 

= a + F5(^) + «- E -r-^"[yiy3]-«-(i + ^)[yiy3] 

n = ^• 

= a + F5( J) - « • [yi ys]-K-[2- AO] • ^ - yf] 

+K ■ [m m] ■ cos [2 • AG] + K ■ ^ - yl] ■ sin[2 ■ A0] 
= a + F5(^) - • [^1 ^3] • {1 - cos[2 • AG]} 

■ I [yl - fi] ■ {sm[2 • AO] - [2 • AO]} . (C.ll) 



C.2 Thin- Lens Transport Map 

Equations (C.9), (C.ll) and (C.3c) finally lead to: 



(C.12a) 



a' — a 



• Ae 



[1 + m)] 



+ 



-n-[x' ■ z']- {l-cos[2- AO]} 



with 



f 

Pa 



+K ■ - [{z'f - {x'f] ■ {sin[2 • AG] - [2 • AG]} ; 
2 



AG = 



K 



Hjso) ■ As 
[1 + fiP^)] 

1 /'(p:) 9 



2[l + /(p^)] 

and Mjq given by (C.IO), where the matrix Aqi appearing in (C.IO) takes the form: 



(C.12b) 
(C.12c) 

(C.12d) 
(C.12e) 



ioi = Ae-H{so) 



/ \ 

-10 



V 0-10/ 



(C.12f) 



56 



(see eqn. (C.7a) ). 

Equation (C.12) contains as special cases the transport maps of a simple solenoid and a 
simple quadrupole already derived in section 4. 
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